Polarization fluctuations in insulators and metals: New and old theories merge 
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The ground-state fluctuation of polarization P is finite in insulators and divergent in metals, 
owing to the SWM sum rule [I. Souza, T. Wilkens, and R. M. Martin, Phys. Rev. B 62, 1666 
(2000)]. This is a virtue of periodic (i.e. transverse) boundary conditions. I show that within 
any other boundary conditions the P fluctuation is finite even in metals, and a generalized sum 
rule applies. The boundary-condition dependence is a pure correlation effect, not present at the 
independent-particle level. In the longitudinal case V • P = — p, and one equivalently addresses 
charge fluctuations: the generalized sum rule reduces then to a well known result of many-body 
theory. 



PACS numbers: 77.22.-d, 71.10.-w, 71.10.Ca 
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In a paper appeared in 2000 Souza, Wilkens, and 
Martin 1] (SWM) proved a fluctuation-dissipation sum 
rule relating the ground-state fluctuation of polarization 
P in a quantum system to its macroscopic conductivity. 
The sum rule implies that P fluctuations are finite 
in insulators and divergent in metals, thus providing 
a clearcut qualitative difference between insulating and 
metallic ground states. In fact SWM complete the 
program initiated in 1964 by W. Kohn with his "Theory 
of the insulating state" 0. However, a different and 
apparently unrelated fluctuation-dissipation sum rule is 
well known since the 1950s in many-body physics 0. |j| 
The disturbing fact is that metals and insulators do 
not behave in a qualitatively different way as far as 
the latter sum rule is concerned. I show here that 
both sum rules are special cases of a more general one, 
the difference owing to the boundary conditions (BCs) 
adopted when taking the thermodynamic limit: SWM 
adopt periodic Born-von-Karman BCs, i.e. transverse, 
while within many-body physics it is customary to adopt 
longitudinal ones. The precise meaning of "transverse" 
and "longitudinal" in the present context is illustrated 
below. For any BC choice different from the purely 
transverse the ground-state P fluctuation is finite even 
in metals. I also show that the BC dependence of 
the P fluctuation is a combined effect of the long 
range of Coulomb interaction and of electron correlation. 
There is no such dependence for independent electrons 
(either Hartree-Fock or Kohn-Sham), where the standard 
SWM sum rule applies anyhow. The presentation starts 
considering a finite system with "open" BCs, and then 
proceeds to taking the thermodynamic limit in the 
appropriate way. 
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Since it is expedient to deal with quantities that become 
intensive in the large- N limit, I define the first and second 
cumulant moments per electron: 
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(Greek subscripts indicate Cartesian components 
throughout). The first moment (times the trivial factor 
—eN/V) is the electronic term in the macroscopic polar- 
ization P, while the second cumulant moment {r a rp) c 
is its quadratic quantum fluctuation in the many-body 
ground state. We notice that the second moment is 
a function of the relative coordinates, while the first 
moment is a function of the absolute ones; indeed, the 
macroscopic polarization P is well defined only when 
the (classical) nuclear contribution is accounted for. As 
said above, SWM address the P fluctuations in extended 
systems within periodic BCs. Therein, the position R is 
a "forbidden" operator [5j and the definition of {r a rp} c 
looks formally quite different from Eq. J3J 6, 7]. 

Indicating with Xj = (r,-,<Xj) the space and spin 
coordinates of the i-the electron, the one-body and two- 
body densities are defined as: 



n(ri) = N 
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rfx 2 • • •dx A r|*(xi,x 2 , . . .xtv)| 2 ; (4) 



Let * be the singlet ground-state wavefunction of an n (2) (ri, r 2 ) = N(N—1) ^ / dx 3 ■ ■ ■ dxjv|*( x i> . . . xat)| 2 
TV-electron system, with even N, within "open" BCs, i.e., CTl(T2 J 

^ is square-integrable. We address the first and second (5) 
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Straightforward manipulations lead to the equivalent 
form: 

{r a rp) c = J drdr' (r-r')a(r-r')/3 

x [n(r)n(r') - n (2) (r, r') ]. (6) 

We observe that, for large values of the relative 
coordinate r — r', the electron distribution becomes 
uncorrelated and one has n^(r, r') ~ n(r)n(r'): 
this fact can be regarded as a manifestation of the 
"nearsightedness" principle |£g. For any finite N the 
integration in Eq. © obviously converges, owing to the 
boundedness of the ground wavefunction 'J. One of the 
main questions is whether (r a r^) c goes to a well defined 
finite limit or instead diverges in the limit of large N. 

A well known exact sum rule relates the two-body 
density to the frequency integral of the imaginary part 
of the linear response: such relationship belongs to 
the general class of fluctuation-dissipation theorems 0, 
Ho| . By definition, the linear polarizability tensor 
a0 7 (w) yields the /? component of the dipole d linearly 
induced by an electric field Eo of unit magnitude in 
the 7 direction, at frequency to. I address purely 
electronic response, therefore assuming clamped nuclei. 
Furthermore I stress that Eo is the field far outside the 
finite sample, different from the screened macroscopic 
field E inside. Starting e.g from Eq. (2.17) in Ref. [n) 
and using Eq. © , it is straightforward to prove the sum 
rule: 
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It is worth noticing that the rhs is by construction a 
ground-state property, while the lhs is a property of 
the excitations of the system. In Eq. J7J|, fluctuation 
and dissipation are perspicuous: from the definition of 
Eq. @ the rhs is a quantum fluctuation, while the 
imaginary part of a(ui) measures dissipation in the zero- 
temperature limit [9j. Because of the sum rule, the actual 
value of (r a rp) c in a given system can in principle be 
measured by actually probing the excited states. 

I now discuss Eq. (|7|) in the limit a macroscopic solid, 
where for the sake of simplicity the bulk is assumed as 
macroscopically homogeneous and isotropic. Therefore 
the macroscopic polarization P = d/V linearly induced 
by a macroscopic field at frequency lo can be written as: 
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where s(u>) is the electronic (clamped-nuclei) macroscopic 
dielectric function of the bulk material. In terms of Eq, 
this polarization is 



P = ^^^T^^T' 



(9) 



where the relationship between E and Eo depends on the 
shape of the sample. 

Following a well-known practice for the study of 
dielectric bodies, we consider a sample of ellipsoidal 
shape, in which case the E field is constant in the bulk, 
and the shape effects are embedded in the depolarization 
coefficients np, with X)/3 n /3 = The main relationship 
is [13: 
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The extremely prolate ellipsoid (n x 



(10) 

= 1/2, n z = 0) 

is a cylinder along z, while the extremely oblate one 
\Tb X - - fly = 0,n z = 1) is a slab normal to z. The 
slab geometry epitomizes both the longitudinal and the 
transverse cases: P is purely longitudinal when along 
z, and purely transverse when along xy. In the former 
case, in fact, we have P z = P z {z) (independent of xy): 
hence V-P^O, V x P = 0. Conversely in the latter 
case we have P x = P x (z) (independent of xy): hence 
V • P = 0, V x P / 0. It is worth noticing that the 
charge is uniquely related to P via V • P = — p in the 
longitudinal case, whereas the charge does not enter a 
macroscopic description in the transverse one. 

In the ellipsoidal geometry Eqs. © and lfTU)l yield [l^ : 
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The a tensor is diagonal over the ellipsoid axes, and the 
above results transform Eq. (J7J) into: 
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It is expedient to recast this sum rule in terms of 
the conductivity o-(lu), which by definition measures the 
macroscopic current linearly induced by a field E at 
frequency u>. Since the current is the time derivative 
of the electronic polarization, Eq. (JHJ yields e(w) — 1 = 
ATtio~{oj) / lj and 



V h 
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(14) 

which generalizes the SWM sum rule. In fact the 
assumption of periodic BCs corresponds — as I am going 
to explain below — to the choice np = 0, yielding the 
original SWM sum rule: 



V h [°° dcu 



(15) 



The rhs has a qualitatively different behavior in insulators 
and in metals. In the latter materials, in fact, the real 
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part of the conductivity is either finite or divergent in 
the dc (lo — * 0) limit, thus implying in both cases the 
divergence of the integral, ergo of the P fluctuation. In 
insulators, instead, the integral in the rhs of Eq. 115|) 
converges to a finite value. 

Some precursor work, before SWM, attempted to 
relate ground-state fluctuations to dc conductivity 0, 
Q] . It is worth noticing that such work severely overlooks 
the role of BCs, while instead SWM provide a rigorous 
theory in a purely transverse framework. However, SWM 
neither consider different BCs, nor relate their work to 
the Nozieres-Pines early fluctuation-dissipation sum 
rule. Here we provide a generalization of SWM to all 
possible BCs, Eq. (fT4"j) . Its novel outstanding message is 
that for any np ^ choice the P fluctuation is finite even 
in metals (contrary to what stated in Ref. 0]). 

The second cumulant moment (r / 3r 7 ) c has been defined 
as a bulk property of the condensed system, which 
measures the quadratic quantum fluctuations of the 
polarization in the many-body ground state at zero 
temperature . It may appear therefore disturbing that 
its expression, as given in Eq. I|14|l . depends explicitly — 
via the np coefficients — on the shape which has been 
chosen for taking the large- 1^ limit. In fact, this is a 
real physical effect and has a simple interpretation. 

The fluctuating polarization P induces a surface charge 
at the boundary of the sample, which in turn generates 
a homogeneous depolarizing field E, which counteracts 
polarization: in the unperturbed Eo = case Eq. I|10|) 
reads 

E p = -AwnpPp. (16) 

Seen in this way, the effect obviously does depend on 
shape. But in condensed matter physics one tries to steer 
clear from any shape issue, and therefore one interprets 
Eq. (|16fl . for any choice of np, as a choice of BCs for 
performing the thermodynamic limit. Indeed, Eq. I|16(l 
becomes the basic one, and any reference to shape is no 
longer needed. The condition J2p np — 1 is not needed 
either. When adopting the usual periodic BCs in all three 
Cartesian coordinates, we are effectively imposing E = 0, 
i.e. n x = n y = n z = 0. From what said above, and 
from Eq. (|14|l . one would obtain the same fluctuation 
when working in a slab geometry and addressing the P 
component parallel to the slab, i.e. transverse. The other 
extreme case of Eq. (|16fl . namely n x = n y = n z = 1, is 
also well known in condensed matter physics. In fact, 
the BCs for zone-center phonon modes in cubic binary 
crystals are E = for transverse modes, and E = — 47rP 
for longitudinal ones [H [H H [3 

There is a complete analogy between the ground-state 
fluctuations of polarization in a many-electron system at 
zero temperature, as discussed here, and the equilibrium 
fluctuations of polarization in a classical dipolar system 
at finite temperature. In the latter case, in fact, it 
is well known [2(| that different BCs lead to different 



fluctuations but to the same value for the static dielectric 
constant, provided the correct fluctuation formula is used 
for each case [2]]]. The shape-dependence is a combined 
effect of interparticle correlations and of the long-range 
nature of the interactions. The analogy goes further, 
since even in the quantum case the dependence on shape 
(or equivalently on BCs) is a pure correlation effect, 
not present at the independent-electron level (cither 
Hartree-Fock or Kohn-Sham), where the many-body 
wavefunction is a Slater determinant. In fact, the second 
cumulant moment (r / gr 7 ) c , Eq. (jSJl, is a function of the 
two-body density: the latter, for the special case of a 
single-determinant wavefunction, is an explicit function 
of the one-body density matrix. As such, it can only 
be affected by the mean E field (i.e. zero, for the 
unperturbed system), and not by its fluctuations. 

This is confirmed by the present sum rule. Starting 
from Eq. (JTJ) , we notice that when we evaluate the rhs us- 
ing the independent-electron two-body density, we must 
interpret the a(uj) tensor in the lhs as the independent- 
electron polarizability, which by construction neglects 
self-consistency effects. Therefore E = Eo i.e., after 
Eq. {TTJ, np = 0. Therefore Eq. {H} reduces to Eq. l(T5|l . 
which is manifestly shape-independent (or BC indepen- 
dent). Incidentally, the conductivity o~(lu) therein must 
be understood as the independent-electron conductivity. 

I now address the special form taken by Eq. I|13|) in the 
longitudinal case, where = 1: the diagonal 

zz component is 

Ti V f°° 1 
( z )c = -TT»m / ^Im— . (17) 

One would obtain the same fluctuation working in a 
slab geometry and addressing the fluctuation of the P 
component normal to the slab. Eq. (|17|l applies to 
correlated wavefunctions, and is invalid for independent- 
electron ones; it provides a finite value both in insulators 
and metals. We are going to verify the above general 
findings on the simplest metal of all, namely, the 
homogeneous electron gas, showing that (z 2 ) c is infinite 
in the noninteracting case, and finite in the interacting 
one. 

In order to make contact with the electron-gas 
literature, we need to introduce the static structure factor 
defined as 

S(k) = i(*|X;e <k -(«-^|*>. (18) 

This is identically expressed in terms of the one- and 
two-body densities as: 

S(k) = 1 + 1- J drj dr' e* k -( r - r 'W 2 )(r,r'), 

= l + ^Rk)| 2 (19) 
+ 4 fdrf ( f'e lk - (r - r ' ) [n^(r,r')-n(r)7i(r')]. 
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For k = the second term is equal to N, which 
obviously diverges in the thermodynamic limit: such 8- 
like singularity is neglected as usual. We then expand 
in powers of k imposing centrosymmctry: therefore the 
second term is quartic, and we have to second order: 

S(k) - - £ ^ / dr f dr> (r - r%(r - r') fj 

a/3 

x [n^(r,r') -n(r)n(r')] 

- '^2( r arp) c k a kp, (20) 

a/3 

Therefore for an isotropic system 

(z 2 ) c = lim S(k)/k 2 . (21) 

For the noninteracting (either Hartree-Foch or Kohn- 
Sham) electron gas the one-body density, and hence S(k), 
are known exactly 0: this in fact leads to a divergent 
Eq. 1)21 [I. Polarization fluctuations are indeed BC- (or 
shape-) independent, and diverge even in the longitudinal 
case, thus confirming our general finding. 

In the interacting case S(k), as defined here, depends 
on shape via Eq. (|20|l . whereas in the existing electron- 
gas literature S(k) is apparently shape- independent. 
The reason is very simple: such literature addresses 
charge fluctuations, not polarization fluctuations. It 
has been stressed above that no macroscopic charge is 
associated to transverse polarization fluctuations: charge 
fluctuations manifest themselves only within longitudinal 
BCs, which are therefore implicitly assumed by electron- 
gas theorists. Our longitudinal Eq. I|17(l . together with 
Eq. yields 

which indeed is the standard fluctuation-dissipation 
theorem for the interacting electron gas, known since the 
1950s 4j. The frequency integral is finite: replacement 
into Eq. i|21[) confirms that the longitudinal polarization 
fluctuation (z 2 ) c is finite as well. 

In conclusion, I have reconciled two different forms of 
the fluctuation-dissipation sum rule for quantum many- 
body systems: one recent and one old jj|. The two 
were apparently contradictory and apparently unrelated. 
Instead, I have shown that a more general sum rule holds, 
yielding the previously known ones as special cases. At 



the root of the generalization is a careful treatment of 
electron correlation in Coulomb systems. Remarkably, 
the novel feature found here is a pure correlation effect, 
not present at the Hartree-Fock or Kohn-Sham level. 
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